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ABSTRACT 

We investigate the optical variability of 7658 quasars from SDSS Stripe 82. 
Taking advantage of a larger sample and relatively more data points for each 
quasar, we estimate variability amplitudes and divide the sample into small bins 
of redshift, rest-frame wavelength, black hole mass, Eddington ratio and bolo- 
metric luminosity respectively, to investigate the relationships between variabil- 
es ity and these parameters. An anti-correlation between variability and rest-frame 
wavelength is found. The variability amplitude of radio-quiet quasars shows 
almost no cosmological evolution, but that of radio-loud ones may weakly anti- 
correlate with redshift. In addition, variability increases as either luminosity or 
Eddington ratio decreases. However, the relationship between variability and 
black hole mass is uncertain; it is negative when the influence of Eddington ratio 
is excluded, but positive when the influence of luminosity is excluded. The intrin- 
sic distribution of variability amplitudes for radio-loud and radio-quiet quasars 
are different. Both radio-loud and radio-quiet quasars exhibit a bluer- when- 
brighter chromatism. Assuming that quasar variability is caused by variations of 
accretion rate, the Shakura-Sunyaev disk model can reproduce the tendencies of 
observed correlations between variability and rest-frame wavelength, luminosity 
as well as Eddington ratio, supporting that changes of accretion rate plays an 
important role in producing the observed optical variability. However, the pre- 
dicted positive correlation between variability and black hole mass seems to be 
inconsistent with the observed negative correlation between them in small bins of 
Eddington ratio, which suggests that other physical mechanisms may still need 
to be considered in modifying the simple accretion disk model. 



Subject headings: galaxies: active — galaxies: nuclei — quasars: general - 
quasars: variability 
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Introduction 



Variability is one of the major characteristics of Active Galactic Nuclei (AGNs), whose 
luminosities vary in all the bands from 7-ray to radio, on timescales from hours to years. Al- 
though the study of variabilit y plays an important role in investigat ing the nature of the com- 
pact central region in AGNs (IHook et al.lll994l ; IWilhite et al.ll2008l ). the mechanism underly- 
ing quasar variability is still i nconclusive. Several physi cal mechanisms have been proposed 
to e xplain quasar variability (Vanden Berk et al.l 2004 ), including accretion disk instabili 



ties 



1992 



mg 



Reeslll984l ; iKawaguchi et al. 



1998 



Kato et al. 



1996 



Manmoto et al. 1996 



2008) , Poissonian processes, such as multiple supernovae or star collisions (jTerlevich et al. 



Czernv et al. 



jHawkinslll99~3h 



Courvoisier et al.lll996l ; iTorricelli-Ciamponi et al.ll2000l ). and gravitational microlens 



To clarify the nature of quasar variability, most previous studies focused on the depen- 
dencies of the variability indicator on redshift, time lag, rest-frame wav elength and lumi- 



nosity A decrease of variability was found towards longer wavelengths ( ICutri et al.l 11985 



Paltani fc Courvoisierl Il994j ; Idi Clemente et al.l Il996l ) , but the i nfluences of other p a rame- 
ters on the correlation were not excluded due to small samples. IVanden Berk et al.l (120041 ) 
confirmed such an anti-correlation by excluding the dependencies of variability on redshift 
and luminosity, however, for the scarcity of data points of each quasar, they estimated one 
variability indicator from a sample of quasars (the ensemble variability method) rather than 
from an individual quasar (the individual variability method). This ignores the diversity of 
an individual quasar light curve, which may eventually fade out the actual relations. Due 
to the statistical incompleteness of small samples, previous studies fou nd controversial ev- 



idences for the dependencies of variability on redshift and luminosity ( Bono 



Netzer fc Shefferlll983l : ICristiani et al.lll990l : iGiallongo et al.lll99ll : IHook et al 



i et al 



1979 



1994J). Thus, 



the results may not be statistically reliable due to possible dependencies of variability on 
other parameters. Only recently, with larger samples, some work managed to disentangle 
the anti-correlation between the variability ampl i tude and luminosi ty by restricting other 
parameters in small ranges (IWilhite et al.l 120081 ; iBauer et al.l 120091 ). However, they still 
ad opted the ensemble variability method to estimate the variability indicator. Moreover, 



as IWilhite et al.l (120081 ) only studied quasars with redshift larger than 1.69, the situation 
for quasars with lower redshifts were still unknown. Besides luminosity, more attention has 
been paid to investigate the influences on variability from other intrinsic parameters, i.e., 
black hole mass and Eddingt on ratio (the ratio of bolometric luminosity to Eddington lu- 
minosity). I Wold et al.l (120071 ) found a correlation between variability and black hole mass 
from about 100 quasars, without accounting for the influences of other parameters. Based 
on t heir ensemble variab ility study of nearly 23,000 quasars from the Palomar QUEST Sur- 



vey, 



Bauer et al.l (120091 ) found the dependencies of quasar variability on black hole mass 
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(positively), and Eddington ratio (negatively). Similar conclusions were reached by lAi et al. 
( 120101 ) based on the individual variability indicators but without disentangling the relations 
between these parameters. 

Although variability depends on many parameters, the underlying physical mechanisms 
are still uncertain. Dividing the sample according to different parameter bins is essential 
to analyze their influences on variability. Moreover, to avoid missing the diversity of the 
individual quasar light curve, the detailed variability indicator of each quasar rather than 
the ensemble variability indicator from a group of quasars, should also be utilized. However, 
no previous work has employed the two essential methods simultaneously in their correlation 
analysis. 

In this paper we revisit the correlations between the variability amplitude and physical 
parameters based on a sample of 7658 quasars in SDSS Stripe 82, with both the individual 
variability method and the detailed parameter binning techniques. We describe the sample 
selection and quasar variability estimations in § [2j The variability distribution and color- 
brightness relationship are given in § [31 § El shows the correlation between the variability 
amplitude and rest-frame wavelength for radio-loud and radio-quiet quasars. § [5] shows the 
correlation between variability and redshift. We disentangle the dependencies of variability 
on Eddington ratio, luminosity and black hole mass in § |6j Implications of our results and 
the comparison with the accretion disk model are presented in § We summarize our results 
in §EJ The cosmological parameters of (ho, Qo, A ) are set to be (0.73, 0.27, 0.73) throughout 
the paper. 



Quasar sample and variability estimation 
2.1. Sample 



A sample of 9254 variab 



e quasars is obtained from cross-matchin g: the spectroscopic con- 



firmed SDSS DR7 quasars ([Schneider et al.ll2010l ; IShen et al.ll201ll ) and a sample of 67507 



variable sources in SDSS Stripe 82, which lies along the celestial equator in the Southern 
Galactic Hemisphere (22h 24m < a J20 oo < 04h 08m, -1.27° < 5, m oo < +1.27°, ~ 290 deg 2 ) 
and have repeated photometric observations (at least 4 per band, with a median of 10) mea- 



sured in up to 10 years in the u'q' r'i'z' system (IFukugita et al.l Il996l ; iMacLeod et al.l 12010 



Ivezic et al. 


2007; 


Sesar et al. 


2007) 



ton ratios are obtained from the quasar catalog in IShen et al.l (120111 ) , where the black hole 



mass es of quasars were settled by fid ucial virial mass est imates: H/3 ( jVestergaard fc Wilkes 



20061 ) estimates for z < 0.7, [Mgll] f lShen fe Kelly! 12010( 1 estimates for 0.7 < z < 1.9 and 
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ClV] ( jVestergaard fc Wilkeal2006l ) estimates for z > 1.9. If more than one e missi on line are 
available, we still adopt the fiducial virial black hole mass in lShen fc Kellyl ( 120101 ). In total, 
the quantity is measurable for 9148 quasars (98.9%) in the sample. 



2.2. Variability estimation 



Ai et al. 


(2010 


) and 


Sesar et al. 


(2007) 
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where n is the number of observations for a single source in each band, m\ is the magnitude 
from the i t h observation, (m) is the mean magnitude, and £ is estimated from the photometric 



error Am;: 



1 " 

- E ( Am ^ 



(3) 



We adopt V calculated in this way as the variability amplitude, which describes the 
intrinsic quasar variability by excluding the photometric error. Before estimating the vari- 
ability amplitude in one band for an individual quasar, we eliminate the photometric out- 
l iers to make sure the calculated variability amplitude is not caused by non-physical reasons 
jSchmidt et all boiol : iMacLeod etaD boiol : butler fc Blooml boilh . Assumi ng that a sig- 
nifica nt decrease in magnitude in a single observation may be non-physical (jSchmidt et al. 
2010l ). we exclude data points with very faint luminosity, i.e., magnitude fainter than 30. 
Data points with large photometric uncertainties (errors larger than 1 magnitude) are also 
excluded. Photometric measurements are marked as outliers if the residuals between them 
and their median magnitude of all the preliminarily cleaned data points are above 5 times 
their standard deviation. This procedure is iterated until all the outliers are excluded, each 
time removing only the strongest outlier and updating the median magnitude along with the 
standard deviation for the remaining measurements. After eliminating all the assumed non- 
physical outliers, if more than 10 cleaned data points remained, we calculate the variability 
amplitude with the equations mentioned above. Otherwise, variability amplitudes are set 
to be 0. To get more reliable statistical results, in the following studies we only investigate 
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quasar variability amplitudes in the SP SS q', r' and i' ba nd, as photometric errors in the 
u' and z' band are substanti ally larger (ISesar et al.l 120071 ). Among the 9148 quasars with 
measurable black hole mass ( Shen fc Kelly 201oh . 8967 quasars have reliable variability in- 
dicators in the g', r' and %' band, which means that variability in each band for a quasar is 
a non-zero value estimated from at least 10 data points, even after eliminating photometric 
outliers. 



3. Variability distribution and color changes 



Since the ph ysical origin of opt ical variability in radio-quiet and radio-loud quasars 
may be different (IGiveon et al.lll999l ). among the 8967 quasars, we divide the 7658 quasars 
with either radio d etections or in the fields of the FIRST radio survey into two subsamples 
( [White et al.lll997l ): the radio-quiet subsample contains quasars with radio loudness below 
10; the radio-loud subsample holds quasars with radio loudness above 10. Here radio loudness 
is defined as the ratio of the ob served radio flux density at rest-frame 6 cm and the optical flux 
density at rest-frame 2500 A( Becker et al.l Il995t IShen et al.l 12011 ). 416 radio-loud sources 
and 7242 radio-quiet sources are obtained for the study in the following sections. 

From Fig. [TJ we can see that the variability amplitudes of the radio-quiet and the radio- 
loud quasars are mostly distributed between 0.05 mag and 0.3 mag, except that there is a 
relatively higher fraction of the radio-loud quasars (~ 2.96%) with variability amplitudes 
larger than 0.3 mag compared to that of the radio-quiet quasars (~ 1.73%). The median 
variability amplitude of the radio-quiet objects in the g', r' and i' band are 0.125±0.071 mag, 
0.106±0.060 mag and 0.095±0.057 mag and that of the radio-loud objects are 0.127±0.096 
mag, 0.106±0.081 mag and 0.096±0.078 mag respectively. Therefore, the median variability 
amplitude of the radio-quiet objects is very similar compared to that of the radio-loud 
quasars. However, their distributions differ significantly according to the results of 100 
Kolmogorov-Smirnov (KS) tests (the median p ~ 0.012). Each time we randomly split 
both the radio-quiet and radio-loud quasar samples into half and perform KS tests on the 
generated samples. For variability values in the g', r' and i' band, 99.7, 91.7 and 81.5 percent 
of p values from these KS tests are less than 0.05 respectively, further indicating significant 
differences between their intrinsic distributions. Moreover, for both these radio-quiet and 
radio-loud quasars, the median variability amplitude decreases systematically as the observed 
wavelengths increase (see §H]for more details). 

Here we take the magnitude difference between two adjacent SDSS filter bands as the 
color indicator. Using the data in the g', r' and %' filter band, we can get 2 colors, namely 
g'-r' and r'-i'. If both the difference of time in which the quasar luminosity of two adjacent 
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bands show their brightest state and the difference of time when luminosity of two adjacent 
bands show their faintest state are less than 2 days, then these quasars are selected to study 
the quasar color behavior. For the color g'-r' and r'-i', we find 2811, 2553 radio-quiet quasars 
and 166, 158 radio-loud ones respectively. The differences between colors in the brightest 
state and faintest state are shown in Fig. |2j Both the radio-quiet and radio-loud quasars 
display a strong bluer-when-brighter chromatism; namely 82.9%, 67.8% of the radio-quiet 
quasars and 77.1%, 63.3% of the radio-loud quasars become bluer in the g'-r' and r'-i' color 
respectively, when they become brighter. 



4. Dependence on the rest-frame wavelength 



As quasar variability amplitudes depend on quasar bolometric luminosity (Lboi), rest- 
frame wavelength (A r f), black hole mass (Mbh) and possibly redshift (z), such dependencies 
can be formulated as V(X r f, z, M B h, -^boi), where the combination of dependencies of the 
variability amplitude on Mbh and Lboi can be replaced by the combination of its depen- 
dencies on Mbh and -Redd (Eddington ratio) or that on -Redd and Lboi- Dependencies of 
the variability amplitudes on each parameter can be obtained in subsamples where other 
parameters are constrained within small ranges. Some subsamples or sub-subsamples ob- 
tained from the combinations of those parameters yield unphysical (imaginary) correlation 
results. In most cases when this occurs the number of quasars is less than 10. In order 
to obtain statistically significant results, we require at least 10 quasars in a data set before 
including it in the analysis. For all further analysis, binned data sets including fewer than 
10 quasars are rejected. The bin size of z is 0.3 and that of -Redd, Mbh and Lboi are all 0.5 
dex. The minimum values of the four parameters are respectively 0.08, 10 70 M Q , 10 44 3 erg/s 
and 10~ 313 , while the maximum values are 5.08, 10 10,8 M Q , 10 48,0 erg/s and 10 0,7 . In this 
work, the ranges of the four parameters are chosen to be [0.08,4.9], [10 75 M Q , 10 10 ' 5 M ], 
[10 44,5 erg/s, 10 48 erg/s] and [10 -2,0 , 10 0,5 ] respectively, as shown in Table [H More than 
99.5 percent of quasars can be covered in this way and further expansion of the ranges does 
not yield more qualified subsamples. Supposing that X is one of the five parameters (A r f, z, 
Mbh, -^boi and -Redd), Spearman's Correlation Coefficient rl_X is calculated in each corre- 
sponding subsample or sub-subsample to describe the rank correlation, while the significance 
p_X of its value deviating away from is also obtained. To parameterize the observed re- 
lation between the variability amplitude and the parameter X, we fit all the data in each 
subsample with a simple function: 



V(X) = bJXlogX + aJC 



(4) 
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Accompanied by the slope 6_X and y axis intercept value a_X, Pearson Product-Moment 
Correlation Coefficient r2_logA" is also calculated to illustrate the strength of the proposed 
relation (r2_X for short). In the following sections, X will be replaced with A r f, z, Mbh, 
Lboi and -Redd- 

To isolate the dependence of the variability amplitude on rest-frame wavelength for the 
radio-quiet quasars, we first divide the 7242 radio-quiet quasars into small bins in a two- 
dimensional space of redshift and Eddington ratio, as shown in Fig. [3j Ranging from the 
minimum value to the maximum value, z is divided into 16 bins, log -Redd is divided into 5 
bins, as listed in Tabled] This procedure yields 16x5 possible subsamples. Each subsample 
is further divided in 6 log Mbh bins. As seen in Fig. [3] 49 subsamples containing more than 
10 quasars in each, which are filled with red points, are further subdivided in log Mbh- The 
indices of the 49 qualified subsamples are shown in the lower corner of the corresponding 
grids. 104 qualified sub-subsamples with more than 10 quasars in each are finally constructed, 
e.g., indices like l_(l-2) means that 2 qualified sub-subsamples are selected from 6 small bins 
after the division of log M B h of the subsample (0.08 < z < 0.4 and— 1.5 < log .Redd < —1.0). 
We calculate rl_A r f, r2_A r f and fit all the data in each qualified sub-subsample with the 
function like Eq. (JIJ). 

To demonstrate the rest-frame wavelength A r f dependence of the variability amplitude V 
with significant statistics, in Fig. H] we show the results of the 25 sub-subsamples containing 
the most numerous quasars, starting with the largest sub-subsample. To make the relation 
more obvious in each panel, A r f of the qualified sub-subsample is evenly divided into 10 
bins. The median variability in each bin is denoted as red diamonds. Error bars show the 
standard deviation of variability values for quasars in each rest-frame wavelength bin. For the 
radio-loud quasars, we utilize the same mechanism to investigate the relationship between 
V and A r f. 8 qualified sub-subsamples are obtained from binning the 416 radio-loud quasars 
in log .Redd, z and log Mbh- The ranges of these binned parameters and the correlations 
between V and A r f are shown in Fig. [21] in the Appendix. Note that for radio-loud quasars, 
figures like Fig. [21] which describe the dependencies of variability on the quasar parameters 
will be shown in the Appendix. For both these radio-quiet and radio-loud quasars, the 
anti-correlation between V and A r f is evident — quasars appear to be less variable at larger 
rest-frame wavelengths. It is further confirmed by the distributions of calculated rl_A r f, 
r2_A r f, b_\ v{ and & err _A r f values , as shown in Fig. [5] In the left panel, the dots on the left 
of rl_A r f = suggest that the majority of the qualified radio-quiet sub-subsamples exhibit 
negative rl_A r f, with a median value of -0.229. The significance p_A r f of its value deviation 
away from is 0.019, which confirms the moderate anti-correlation at a confidence level over 
98%, while the distribution of the asterisk signs presents that the median rl_A r f value for all 
the 8 radio-loud sub-subsamples is -0.187 with a lower significance level (p~0.218). In the 
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middle panel, r2_A r f values for the 100 radio-quiet and all the 8 radio- loud sub-subsamples 
are smaller than 0. The median values are -0.233 and -0.224 for the radio-quiet and the 
radio-loud data sets respectively, while the distribution of the radio-quiet quasars is broader 
than the radio-loud data sets. From the fitting procedure, the median values of 6_A r f and 
a_A r f for the radio-quiet sub-subsamples are -0.137±0.052 and 0.579±0.180, while the median 
6_A r f and a_A r f value for the radio-loud sub-subsamples are -0.225±0.143 and 0.950±0.494. 
The relative ratios of 6_A r f and their errors 6 orr _A r f are shown in the right panel, where the 
ratios are generally larger for these radio-quiet sub-subsamples, implying a more robust linear 
relation. 



5. Dependence on redshift 

We first exclude the influences of the intrinsic quasar parameters on the variability 
amplitude by dividing the radio-quiet and radio-loud quasar sample into small bins in the 
log MeH-log -Redd space with their bin sizes of 0.5 dex. For radio-quiet quasars, 23 subsam- 
ples are considered qualified, with more than 10 quasars in each one. As shown in Fig. [6j 
23 grids marked with algebraic numbers in the lower side of the grids and filled with red 
points refer to those qualified subsamples. For the radio- loud quasars, only 11 subsamples 
are qualified. To eliminate the influence of rest-frame wavelength on the relation between 
the variability amplitude and redshift, we prefer to divide A r f into small bins. However, as 
the rest-frame wavelength and redshift are closely related, the binning of A r f in one indi- 
vidual SDSS filter band will lead to a decrease of the redshift range, causi ng the evaluated 



dependence of variability on redshift unreliable ( iVanden Berk et al.l 120041 ) . Therefore, we 
treat variability amplitudes in the 3 bands (</, r' and i') for an individual quasar as vari- 
ability amplitudes for 3 individual quasars with the same redshift. Based on the variability 
amplitude and Arf of the triple-enlarged subsample, we choose an intersecting range of Arf 
in the 3 bands for quasars in each qualified subsample, which are binned with an equal 
size of 400 A. In this way, we manage to exclude the influence of A r f on the correlation 
between V and z by restricting the sub-subsample in small ranges of rest-frame wavelength, 
but without decreasing the range of z. We use the same method mentioned in § H]to select 
qualified sub-subsamples, e.g., indices like l-(l-4) in Fig. [6] means that after the A r f division 
of the subsample (10 7 - 5 M < M B h < 10 8 M Q and -1 < log .Redd < -0.5) into several 
sub-subsamples, 4 of them are qualified with more than 10 quasars in each. The ranges of 
M B h and -Redd for each sub-subsample are shown in Fig. E] and the range of A r f is shown 
further in the lower right of each panel in Fig. [71 



The relation between the variability amplitude and redshift for the 25 qualified radio- 
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quiet sub-subsamples containing the largest numbers of quasars are shown in Fig. including 
the values of rl_z, r2_z, b_z and the ranges for all the restricted parameters. We evenly 
separate the quasars into 10 redshift bins, to make these tendencies more obviously seen. 
The median variability amplitude V and redshift z in each bin are denoted as red diamonds, 
and error bars show the dispersion of variability values for these quasars. These results 
are fairly noisy, and it is difficult to detect any clear trend with redshift, which is further 
confirmed in Fig. [HJ Again in a similar way, the 11 qualified radio-loud subsamples are 
divided into sub-subsamples with a bin size of rest-frame wavelength of 400 A. 18 qualified 
sub-subsamples are obtained and details of their redshift dependence of variability are shown 
in Fig. [22] in the Appendix. 

rl_z, p-Z, r2_z and the results from the linear fitting for both the radio-quiet and radio- 
loud data sets are shown in Fig. [SJ 78.3% of the radio-quiet sub-subsamples show negative 
rl_z and r2_z values. However, as 42.6% of the negative rl_z values for the radio-quiet ones 
are close to zero within a range of 0.1, the median rl_z value is -0.081 above 81.1% confidence 
level, suggesting almost no relation. Meanwhile, the middle panel indicates that no linear 
relation between them exists, consistent with the distribution of |6_2;|/|6err-^| in the right 
panel — 65.0% of these data sets show |&_z|/|& crr _z| values less than 2. For the radio-loud 
data sets, 61.1 percent show negative relations; the median rl_z value is -0.153 with the 
median confidence level ~ 0.306). The very weak negative relation has to be concerned 
with a larger radio-loud data set due to its low confidence level. Distributions of r2 and 
|6_2?|/|£>err-^| show no linear relation between them — 55.6 percent of these sub-subsamples 
exhibit an anti-correlation between V and z, while the other ones show a positive relation. 



6. Dependence on black hole mass, luminosity and Eddington ratio 

Both the radio-quiet and radio-loud sample are first divided into 16 redshift bins as 
mentioned in §|H For the 14 qualified subsamples of radio- quiet quasars with redshift ranging 
from 0.4 to 4.3 and 9 subsamples of the radio-loud quasars with redshift spanning from 
0.4 to 3.1, we calculate Spearman's Rank Correlation Coefficients between the variability 
amplitude in the 3 SDSS bands (g', r' and i') and the three intrinsic quasar parameters 
(Mbh, -^boi and -Redd), without accounting for the relationships among the intrinsic quasar 
parameters themselves. Results for the radio-quiet and radio-loud quasar subsamples are 
shown as magenta diamonds and triangles in the left panel of Fig. (THJ Fig. [TH] and Fig. [T9"j 
respectively. 

For both the radio-loud and radio-quiet subsamples, there is a weak anti-correlation 
between variability and -Redd at a high confidence level. To make a distinction between 
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the Spearman's Rank Correlation Coefficients here and rl calculated after considering rela- 
tionships among the intrinsic quasar parameters, we denote results here as the format like 
rl n _X. rl n _i?EDD of 12 radio-quiet and all the radio-loud subsamples are less than and 
the median values for these radio-quiet and radio-loud subsamples are approximately -0.170 
and -0.390 respectively. Most of the p"_.Redd values are around 0. A weak anti-correlation 
between the variability amplitude and L bol shows up as more than 90% subsamples exhibit 
negative rl n _L^ \ values at a confidence level over 99%. The median rl n _Lboi values for the 
radio-quiet and the radio-loud subsamples are -0.290 and -0.350 respectively. However, no 
conclusive results are found for the relation with Mbh- This is consistent with the anti- 
correlation between variability and L^ Q \ found in previous studies on smaller samples, where 
the relations among; the intrinsic quasar parameters were not taken into co nsideration either 



( iCristiani et al.lll996l ; Idi Clemente et al.lll996l : iPaltani fc Courvoisierl 119971 ). The correlation 



between variability and Mbh for the 0.08 < z < 0.7 radio-quiet sample (rl n _M-Qu ~ 0.26, 
p n -M BU ~ 0.01 for the quasars within 0.08 < z < 0.4 and rl n _M BH ~ 0.13, p n _M BH ~ 0.01 
for the quasars within 0.4 < z < 0.7), corresponds well t o the signifi c ant d ependence of 



variability on Mbh for the < z < 0.75 sample reported by I Wold et al.l ( 120071 ) . We discuss 
this further in § [71 

However, the relations among the intrinsic quasar parameters have to be taken into 
account to investigate the influences of different parameters on the variability amplitude. As 
more luminous quasars tend to posses more massive black holes, the relationship between 
variability and black hole mass will also express itself in the correlation between variability 
and luminosity or the relation with Eddington ratio. As shown in this section, after consid- 
ering the relationships among the intrinsic quasar parameters, trends are still clear in some 
cases but become ambiguous or fade out in other cases. The radio-quiet quasars and radio- 
loud quasars are subdivided into small data sets by binning redshift and one of the three 
quasar parameters to exclude their influences on the relationships between the variability 
amplitude and the other two quasar parameters. The separate consideration of variability 
amplitudes in the 3 filter bands manage to restrict the rest-frame wavelength of each data 
set to small ranges. 

Fig. [3] shows the details of binning in the log -Redd-^ space for the radio-quiet quasars 
and yields 49 qualified subsamples, as described in § HI Fig. [H] and Fig. [TU] show the details 
of binning in the log M-qh-z space and logLboi-^ space, respectively. In Fig. [9], among 16x6 
grids, 52 grids refer to the qualified subsamples. In Fig. [TOj 16x7 grids are obtained and 
45 grids are qualified. For the radio-loud quasars, 11, 11 and 13 qualified subsamples are 
extracted in the log R-^dd-z, log M B n-z and logL bol -z spaces, respectively. In each qualified 
subsample, we further calculate rl, r2, b and a for the correlations between V and the 
other two unrestricted physical parameters. Their distributions are shown statistically in 
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Fig. H31 [16] and [19j For the subsamples obtained from binning in the log Redd-z space, the 
corresponding signs (plus signs for results from the radio-quiet data sets and asterisk signs 
for results from the radio-loud data sets) are colored as red, while for the subsamples from 
binning in the log M B}i -z and log L ho \-z space, signs are colored as blue and grey respectively. 
Thus in the panels of Fig. [131 Fig- (HI and Fig. [191 which contain plus and asterisk signs, 
the number of red, blue and grey plus signs are 49x3 (the g', r' and %' filter band), 52x3 
and 45x3, while the number of red, blue and grey asterisk signs are 11x3, 11x3 and 13x3 
respectively. 

To demonstrate the dependence of the variability amplitude on Lboi in the qualified 
subsamples obtained from the log -Redd-^ space and log M^u-z space, we further show the 
detailed relations for the 25 largest radio-quiet quasar data sets from the former space in 
Fig. [TT]and that from the latter space in Fig. [T2J To make these underlying tendencies more 
intuitively clear, in each subsample the unrestricted parameter Lboi on which we are studying 
the dependence of the variability amplitude, is further subdivided into 10 bins with an equal 
number of quasars in each one. The median variability amplitudes and median parameter 
values are shown as colored diamonds. The error bars represent the 68 percent dispersion 
of variability amplitudes for the quasars in each bin. Only in the r' band the variability 
amplitudes for all the quasars are further shown as grey points. The fitting of the linear 
relation between V and logLboi for all the data in each band are carried out, but only the 
resulting slope in the r' band is shown in the upper right of each panel. rl_Lb i, P-Lbo\ an d 
r2_Lboi for the quasars in each subsample in the r' band are also listed. 94.0 percent of the 
two figures presents an evident anti-correlation between V and logLboi, consistent with the 
statistical study shown in Fig. [131 

As shown in Fig. [131 almost all subsamples exhibit an anti-correlation between the 
variability amplitude and Lb i after considering the relations among the intrinsic quasar 
parameters. The median value of rl_Lboi for the radio-quiet subsamples within small ranges 
of z and log -Redd is -0.269 over 98% confidence level, while the median value of rl_L bo i 
for the radio-quiet subsamples within small ranges of z and logMeH is -0.343 above 99% 
confidence level. The distribution of r2_Lboi with a median value of -0.260 and -0.335 for 
the subsamples in the log Redd~z space and logM B H-^ space, respectively, as displayed in 
the middle panel, indicates a weak linear relation between V and logLboi- The distribution 
of the fitting slope 6_Lb i and its error is shown in the right panel. Based on the averaged 
value, the anti-correlation between the variability amplitude and Lboi can be described with 
a format similar to Eq. (T4j). The median values of fr_Lb i and a_Lboi for the radio-quiet 
subsamples in the log R^edd - ^ space are -0.044±0.020 and 2.128±0.913, and for the radio- 
quiet subsamples in the logM B H-^ space are -0.063±0.020 and 3.032±0.926. While for 
the radio-loud subsamples in the log Ret>t>-z space, the median rl_Lb i, ?"2_Lboi, &-Lboi and 
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a_Lboi are -0.266 with a lower confidence level (p_L bol ~0.297), -0.269, -0.058±0.054 and 
2.794±2.456; and in the log Mbu-z space, the median rl_L bol , r2_L bol , b_L ho \ and a_L bo i are 
-0.384 with a significance of 0.164, -0.321, -0.071±0.055 and 3.373±2.536. As displayed in the 
right panel, the error in b_L hol is comparable to b_L ho y for the radio-loud quasars, indicating 
that the dispersion is large and the linear relation is very weak. According to the simplest 
Poissonian model, in which starburst events are independent and dis crete, the luminosity 



varia bility is expected to vary with luminosity as logV oc 0.5 logL bol ( iCid Fernandes et al. 



20001 ) . Despite the possibility that a wide range of b values could be obtained by considering 



more detailed Poisso nian models, where sup ernovae and their remnants were considered, b 



is oft en close to -0.5 fjVanden Berk et al.ll2004t iPaltani fc Courvoisierl 119971 ; lAretxaga et al. 



19971 ). We fit our data with a two-parameter function log V = b log L bo i + a and the median b 
is -0.223±0.075. Thus, the results obtained from our sample does not support the Poissonian 
model. 

The dependence of the variability amplitude on -Redd is analyzed based on the 45 
subsamples in the logL bo i-z space and 52 qualified subsamples in the logMeH-^ space. To 
specify the relation, the details for the 25 largest subsamples in each space are further shown 
in Fig. [T3] and Fig. [T5] respectively. Going from the upper to the lower panels within various 
parameter ranges in Fig. [T5], the anti-correlation between V and -Redd is apparent. While in 
only 4 panels of Fig. [T4J we find no discernible relation or an insignificant correlation, such 
as the third panel in the lowest line. 

Following the analysis of the relationship between variability and luminosity, statistical 
results of rl_R E DD, P-Redd, ^2_i? E DD, b and |6_-REDD|/|fren— ReddI are shown in Fig. [16j 
For the radio-quiet quasars in the logL bol -2; space and log M B}1 -z space, more than 80.0 
percent and 94.2 percent of rl_i?EDD are negative respectively, indicating an anti-correlation 
between V and -Redd- All the qualified radio- loud subsamples from the two spaces exhibit 
negative v 1 _Redd values. Respectively in the two spaces, the median r 1_-Redd value is 
-0.167 (p_Redd ~ 0.137) and -0.303 (p_.Redd ~ 0.013) for the radio-quiet quasars; and - 
0.304 and -0.369 for the radio-loud quasars over a lower confidence level (jo_-Redd ~ 0.300 
and 0.188). As shown in the last 2 panels, 79.3% of the qualified radio-quiet subsamples 
obtained from the log L^ i-z space exhibit negative t2_.Redd and 6_.Redd values, confirming 
the weak anti-correlation between V and log -Redd- The linear relation for the qualified 
radio-quiet subsamples binned from the logMBH-z space is more evident with a median 
r2__R E DD of -0.299. More than 92.3 percent of radio-loud subsamples from the logL bol -2 
space and all the qualified ones from the logMeH-^ space display negative linear relations 
between V and log -Redd, with median t2_.Redd values of -0.274 and -0.291 in the logL bo i-,2 
and log Mbh-^ space respectively. The linear relation is also described with a simple function 
like Eq. (jlj). |6_i?EDD|/|^err--REDD| values are generally larger for the radio-quiet subsamples, 
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indicating a stronger linear relation. The median values of 5_.Redd and gl-Redd for the 
radio-quiet subsamples in the log Lboi- 2 space are -0.020±0.014 and 0.091±0.011, while for 
the radio-loud quasars the results become -0.045± 0.032 and 0.079±0.035 respectively. For 
the radio-quiet subsamples in the log M-qk-z space the median values of 6_Redd and <2_Redd 
are -0.052±0.019 and 0.075±0.013; for the radio-loud subsamples they are -0.069±0.047 and 
0.056±0.055, respectively. 

The relation between the variability amplitude and Mbh for the 25 largest qualified 
data sets obtained from binning in the log -Redd - ^ space is shown in Fig. [TTJ and that for the 
subsamples with restricted Lboi values is shown in Fig. [181 Almost all the subsamples in the 
former figure report an apparent inverse dependence of V on log Mbh- Proceeding from the 
first panel to the last one in the latter figure, a positive relation is shown in general. Most 
subsamples show a significant positive log M B h dependence, especially for those containing 
more quasars (the first 8 panels). The opposite trend shown in some cases (such as the 3rd 
panel in the lowest line) is very weak and there are too few cases to claim deviations from 
the general trend. 

The calculated rl_MBH, P-Mbh, r2_MBH and 6_Mbh values for each subsample are shown 
in Fig. [T9j In the subsamples obtained from binning in the log R^bd-z space, we can find 
that for both the radio-quiet and the radio-loud quasars, there is a significant anti-correlation 
between them with the median coefficient rl_MBH about -0.203 above 94% confidence level 
and -0.220 at a lower confidence level (p_Mbh ~ 0.432), respectively. While for the radio- 
quiet subsamples with Lboi limited, more than 71.1 percent show positive t1_Mbh values 
and its median value is 0.113 (p_M B H ~ 0.235); the median value of 0.246 (p_M B H ~ 0.386) 
is shown for the radio-loud quasars, suggesting an insignificant correlation. Based on the 
median values of the linear fitting results (6_Mbh and glMbh), we can parameterize the 
relation between V and Mbh using Eq. (@J with X as log Mbh, where 6_Mbh = — 0.033±0.019 
and glMbh = 0.394 ± 0.168 for the subsamples binned from the log -Redd- z space; 6_Mbh = 
0.013 ± 0.014 and o_Mbh = —0.002 ± 0.119 for the subsamples given by binning in the 
logLboi-£ space. The median values of &_Mbh and a_M B H are -0.052±0.052, 0.600±0.477 in 
the log -Redd- z space and 0.040±0.035, -0.206±0.309 in the logLboi-^ space for the radio- 
loud quasars. It guarantees the information provided by the distribution of tT_Mbh — the 
relations are insignificant for the radio-loud subsamples. The median rl, p, r2, b, b eTT , 
a and a crr values for the relationships between V and different quasar parameters are all 
summarized in Table [2j 
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Discussion 



7.1. Comparison with previous work 



Vanden Berk et all ( 120041 ) argued that variability of radio-loud quasars are larger than 
radio-quiet ones, while no evidence for a chan ge in the average op tical variability amplitude 
with increasing radio luminosity is shown in iBauer et al. 6x3). Comparing the median 
variability amplitude, we do not find larger variability amplitudes for radio-loud quasars than 
radio-quiet ones. However, their intrinsic distributions are different based on the results of 
KS test. The variability of radio-loud quasars may be caused by some m echanisms other 
than the accretion disk instability, such as the outflow (IGiveon et al.lll999l ). 



We confirm the anti-correlation between the variability amplitude an d rest-frame wave 



length which is mostly studied in the form of ensemble structure functions (ICutri et al.lll985 



Paltani fc Courvoisierlll994t IVanden Berk et al.l I2004T ) . We show that both the radio-quiet 



quasars and radio-loud qua sars exhibit a 



previous studies ( 


Wil 


lite et a 


200 


5 


Schmidt et al. 


2012: 


Ai et al. 


2011 


)• 



jluer - when-brighter chromatism, in good accor- 



Meusineer et al. 


2010; 


Trevese & Vaenetti 


Honma et al. 


( 


1991 


) showed that when in- 



stabilities occur in the inner disk, the integrated color will change because of the outward 



propa gating hot wave, which can account for the bluer-when-brighter behavior. ISakata et al. 



( 120 111 ) also suggested that it may be caused by changes of accretion rate. Besides this expla- 
nation, radiations from diffe rent regions of accretion d i sk or influence s of host galaxy may 
also induce such a behavior (jShieldslll978t lMalkanlll983t lHawkinslll993l ). However, as quasar 
optical luminosities are dominated by their disks rather than their host galaxies, accretion 
disk instabilities are more likely than influences of host galaxies to explain the bluer-when- 
brighter chromatism. Trevese & Vagnetti (2001, 2002) analyzed the dependence of spectral 
changes on the brightness of quasars and suggested that the chromatism is at least consistent 
with the temperature variation of a single emitting blackbody and can be explained by bright 
spots on the d i sk pos sibly produced by the accretion disk instability, which also agrees with 
Schmidt et all ( 120121 ). It is also possible that a combination of these physical mechanisms 
produces the color chromatism. 

We find no relation between the variability amplitude and redshift for radio-quiet 
quasa rs, which seems contrary to the negat i ve or positive corre lations reported in previous 
w ork (ICristiani et al.lll990l ; iHook et al.lll994t iTrevese et al.lll994l ) . The anti-correlation found 
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Crist iani et al 



(119901 ) may reflect the anti-co rrelation between V a nd Lboi, as their sample 
is too small to disentangle the relation. Only ITrevese et al.l (119941 ) binned the magnitude 
limited sample within small luminosity and redshift ranges and found an anti-correlation, 
which is ascribed to the dependence of variability on wavelength. Some previous work 
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( IGiveon et al.l Il999t ISchmidt et all |2012| ) also found no relation between V and z. After 
corre cting the contributio n from emission lines, the color variability is indepe ndent of red- 
shift (ISchmidt et al.l |2 0121 ). which agrees with our results. However, although iGiveon et al. 
( 119991 ) also found no relation between V and z, their conclusion is probably not reliable due 
to the small size and the narrow redshift range of the quasar sample. Our results do isolate 
the influences of rest-frame wavelength, black hole mass and Eddington ratio and suggest 
that the redshift dependence of radio-quiet quasar variability may not exist. The possible 
anti-correlation between V and z for radio-loud quasars is very weak and not significant; it 
may need to be investigated with a larger sample in the future. 

By investigating the relationships between the variability amplitude and the intrinsic 
quasar parameters in different redshift ranges, both before and after binning in one of the 
parameters, we find some results that are different from those in previous wo rk. In spite 
of oby ious selection effects of the flux-limited sample within < z < 0.75, I Wold et al. 
( 120071 1 reported a significant dependence of variability on black hole mass, approximately 0.2 
magnitude increase of the variability value over 2-3 orders of magnitude in black hole mass. 
They used the maximum magnitude difference of light curves as the variability indicator, 
which is approximately 4 times larger than our estimated variability value. Before considering 
the relations among the intrinsic quasar parameters, the quasars within 0.08 < z < 0.4 and 
0.4 < z < 0.7 in our sample also show a relation between the quasar variability amplitude and 
black hole mass. However, this is not the whole story for the correlation between variability 
and Mbh- We do not find any certain relation between the variability amplitudes and black 
hole mass in the subsamples binned in the intrinsic quasar parameters. Most subsamples 
in narrow ranges of Eddington ratio exhibit a weak anti-correlation while the subsamples in 
narrow luminosity ranges show no relation or a weak correlation between variability values 
a nd black hole mass. We derived the best-fit parameters with the function similar to Eq. (10) 
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Bauer et al.l (120091 ) for our subsamples in small bins of luminosity. Although the resulting 



slope for the largest subsample is 0.067±0.018, from statis tical results of all q ualified radio- 
quiet samples we do not find the weak correlation stated in lBauer et al.l ( 120091 ). The median 
slope is 0.051±0.056 compared to 0.13±0.01 in their work. 



Wilhite et al.l (120081 ) also found a correlation between the variability amplitude and black 



hole mass and an anti-correlation between the variability amplitude and luminosity indepen- 
dent of black hole mass. The black hole mass range of our sample (10 7 M Q - 10 10,8 M Q ) is 
larger than that of their sample, which makes our results statistically more significant. Com- 
pared to their results, the relation between the variability amplitude and luminosity is weaker 
in our case. This may be due to the fact that in our sample, many bright (faint) quasars 
also show large (small) variability, which further shallows the anti-correl ation. Furth e rmore , 
this corresponds well to the anti-correlation found in small Mbh bins in lBauer et al.l ( 120091 ). 
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where the logarithm value of the variability amplitude is proportional to (-0.205 ± 0.002) 
logLboi- With our sample, we calculate the slope using the fitting function similar to theirs 
and th e slope value is -0.263=b0.074. We co nfirm the anti-correlation between variability and 
-Redd (IWilhite et al.ll2008l ; lAi et al.ll2010l ) by excluding the influences of other parameters 
on the variability amplitude. 



7.2. Implications of the standard accretion disk model 

The optical continuum of low -redshift quasars can often be well fitte d by the standar d 
accretion-disk model proposed by IShakura fc Sunyaevl fll973[ ). Following iLi fc Cad ( 120081 ). 
we investigate the dependencies of the variability amplitude on rest-frame wavelength and 
the intrinsic quasar parameters, assuming that quasar variability is caused by a change 
of accretion rate, which is the result of accretion disk instabilities. However, we do not 
utilize the temperature distribution shown in their work, as it is actually an approximation 
to the equatorial temperature in the outer region of the disk and depends on viscosity. 
The temperature distribution we used here describes the radial dependence of the effective 
temperature derived from the emergent flux from an accretion disk around a Schwarzchild 
black hole, which is independent of viscosity, namely, 

1 1/4 

T(R) 



8naR 3 



1 



R m /R 



(5) 



where _R r 



3-/?g, Rg 



2GM 



Here, a is the Stefan-Boltzmann constant, R is the radius of 



the disk, M is the black hole mass, M is the accretion rate. After calculating the emitting 
spectrum at each radius, we integrate it over the whole disk to get the final spectrum. The 
emitted spectrum of the disk can be written as: 

47r/icos(/3)i/ 3 RdR 



fu 



R-m 



e hu/kT(R) 



(6) 



where R out is the outer radius of the disk, h is Plank's constant, (3 is the inclination of the 
disk with r espect to the line o f sight and D is the luminosity distance from the observer to 
the quasar (IFrank et al.ll2002l ). We simply assume the inclination to be 0. It is a reasonable 
assumption since any inclination other than will just reduce the normalization of the whole 
spectra, which does not affect the variability amplitude. In addition, the inclination is not 
likely to change significantly. 

To better compare theoretical results with observations, we obtain r\fa values in the 
model according to the representative -Redd values in our sample. 



M = mM CTit = rjM EDD m, 



(7) 
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where M crit = Ledd/c 2 is the critical accretion ra te, and Mppn = Lvr>r> I ivc 2 ) is the Edding- 



ton accretion rate. Based on this equation, fn in lLi fc Caol (120081 ) is 77 ( ~ 1/16) times our 
fn as their fn is the ratio of M to Medd- 

We calculate spectra for different black hole masses and accretio n rates. Mbh equals to 



10 8 , 10 8,7 and 10 9 M with an observational uncertainty of 0.4 dex (jVestergaard fc Wilkes 



200 6]), while -Redd eq uals to 0.0316, 0.1259 and 0.3162 with an uncertainty of 0.3 dex 



(iKollmeier et al.ll2006l ). After convolving with the SDSS filter response function, we con- 
vert the flux to the SDSS AB magnitude ^(Mbh? fn). Assuming the change of fn is xfn, we 
adopt V = A(MBu,m) — A(Mbr, (1 + x)fn) as variability indicators. Results for the three 
x values (0.4, 0.2, 0.1) are shown as red, blue and black lines respectively in Fig. [20j At a 
specific x, the variability amplitudes along with their errors are calculated as below. A set of 
100 normally-distributed random Mbh around 10 8,7 M and a set of 100 normally-distributed 
random -Redd around 0.1259 are carried out within their observational uncertainties. Based 
on the set of Mbh and the set of -Redd, we estimate variability amplitudes in each filter band 
to see its dependence on rest-frame wavelength. The mean value of the set of estimations 
and 3 times the standard deviation in each filter (u', g', r', %' and z') are adopted as the final 
variability and its uncertainty respectively, as shown in the upper left panel of Fig. [2DJ The 
variability amplitude declines with the increasing of rest-frame wavelength under the three 
changes of accretion rate. In the upper right panel, an anti-correlation between V in the r' 
band and -Redd is displayed, where V and the error in V are respectively the mean value and 
3 times the standard deviation of a set of variability amplitudes estimated from 100 random 
realizations of the models. As luminosity is only related to -Redd if is fixed, an inverse 
relation with Lboi is naturally expected. We do see this trend in the lower right panel of 
Fig. [20l For a set of 100 normally-distributed random -Redd values around 0.1259 within 
its observational uncertainty, we obtained 3 sets of V values, each set with respect to one of 
the three Mbh values. The mean value and 3 times the dispersion of V values in each set is 
adopted as the final V and error in V for the Mbh value. As displayed in the lower left panel, 
there is an increase of optical variability with the increase of black hole mass. Comparisons 
between our theoretic results and observed relationships are carried out in each panel. We 
find that the variability amplitude caused by a change of accretion rate around 0.2m agrees 
well with the majority of variability amplitudes in observations. Moreover, besides the re- 
lation of variability with rest-frame wavelength, the model also reproduces the tendencies 
of correlations between the variability amplitude and luminosity, as well as accretion rate. 



Mbh in the model is generated from a normal distribution with a mean value of 10 8 - 7 M© and an 
observational uncertainty of 0.4 dex, rh is one of the three different rh and xrh refers to a change of the 
corresponding rh. 



- 18 - 



We fit the simulated data points as we do with the real data in the x = 0.2 case, which 
are shown in the upper right of each panel. All the anti-correlations are flatter than the 
results from observations. The predicted correlation seems to be inconsistent with the weak 
anti-correlation between the variability amplitude and Mbh for the subsamples within small 
ranges of -Redd as we stated in § [6j It suggests that other physical mechanisms may still 
need to be considered in modifying the simple accretion disk model. 

8. Conclusion 

A sample of 7658 quasars drawn from SDSS-DR7 is used to investigate the relation- 
ship between the variability amplitude and rest-frame wavelength, redshift and the intrinsic 
quasar parameters, i.e., black hole mass, bolometric luminosity and Eddington ratio. We 
analyzed the radio-quiet and radio-loud quasar variability properties with the sample. Be- 
cause we have better quality data points for each quasar and a large sample, the variability 
indicator is calculated for each quasar from its long term light curve and the dependence of 
different parameters are disentangled by binning the parameter space into small slices. Our 
main conclusions are the following: 

1. There is an anti-correlation between the variability amplitude and rest-frame wave- 
length after excluding the influences of redshift and the intrinsic parameters on variability 
amplitude. The quasar variability amplitude of radio-quiet quasars seems not dependent on 
redshift, while the possible anti-correlation between that of radio- loud quasars and redshift 
still needs be confirmed with a larger sample. 

2. Although the median variability amplitude of the radio- loud quasars is similar com- 
pared to that of radio-quiet quasars, their intrinsic distributions are different. It implies 
that there may be some mechanisms other than the accretion disk instability behind the 
variability of radio-loud quasars. 

3. There is a bluer- when-brighter color chromatism for both radio-quiet quasars and 
radio-loud ones. Some mechanisms, including instabilities occurring in the inner disk, the 
influence of the host galaxy, bright spots or a combination of them, may explain the chro- 
matism. 

4. The variability amplitude is significantly anti-correlated with -Redd and Lt, i in 
the subsamples without binning the quasar parameters. There is no certain correlation 
between the variability amplitude and Mbh in all these subsamples. After considering the 
relationships among the intrinsic quasar parameters, we find that small data sets still exhibit 
an anti-correlation between the variability amplitude and -Redd ; an d between the variability 
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amplitude and Lboi with a lower significance level. Moreover, an anti-correlation of the 
variability amplitude with Mbh begins to emerge in bins of -Redd ; but a positive correlation 
with Mbh is shown in bins of Lb i- 

5. Based on the Shakura-Sunyaev model, assuming the change of accretion rate as 
the origin of quasar variability, we reproduce the relationship of the variability amplitude 
with rest-frame wavelength with good agreement, and the tendencies of correlations between 
variability and -Redd as we U a s A,oi- The predicted correlation between the variability 
amplitude and Mbh seems not to correspond well to the case inferred from observations, 
namely an anti-correlation with Mbh in small bins of -Redd- It implies that the change of 
accretion rate is important for producing the observed optical variability but other physical 
mechanisms still need to be considered in modifying the simple accretion disk model for 
quasars. 
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Table 1. The bin details of each parameter 



Parameter 


Minimum 


Maximum 


Bin size 


NbinB 


(1) 


(2) 


(3) 


(4) 


(4) 


z 


0.08 


1.9 


0.3 


16 


log A/ BH (M q ) 


7.5 


10.5 


0.5 


6 


log^boi (crg/s) 


44.5 


48 


0.5 


7 


log Redd 


-2.0 


0.5 


0.5 


5 



Note. — Col. (1) Parameter names. Col. (2) Minimum 
value of the binned parameter in Col. (1). Col. (3) Max- 
imum value of the binned parameter in Col. (1). Col. (4) 
Bin interval of the binned parameter in Col. (1). Note the 
first bin size of z is 0.32. Col. (5) Number of bins. 



Table 2. Correlation between variability and each parameter 



Radio loudness 


Parameter 


Restricted Parameters 


Spearman relation 


Pearson relation 


Linear fit 






xl 


x2 


x3 


rl 


P 


r2 


6 


a 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


radio-quiet 


log A r f 


Redd 


z 


A/bh 


-0.229 


0.019 


-0.233 


-0.137± 0.052 


0.579 ± 0.180 


radio-quiet 


logz 


Redd 


M BH 


A rf 


-0.081 


0.189 


-0.090 


-0.049± 0.044 


0.131 ± 0.009 


radio-quiet 


log L b oi 


z 


Redd 


A r f 


-0.269 


0.015 


-0.260 


-0.044± 0.020 


2.128 ± 0.913 


radio-quiet 


log Lbol 


z 


M bh 


A r f 


-0.343 


0.003 


-0.335 


-0.063± 0.020 


3.032 ± 0.926 


radio-quiet 


log /? EDD 


z 


/-bol 


A r f 


-0.167 


0.137 


-0.147 


-0.020± 0.014 


0.091 ± 0.011 


radio-quiet 


log Redd 


z 


M BH 


A r f 


-0.303 


0.013 


-0.299 


-0.052± 0.019 


0.075 ± 0.013 


radio-quiet 


log A/bh 


z 


^bol 




0.113 


0.235 


0.104 


0.013± 0.014 


-0.002 ± 0.119 


radio-quiet 


log A/bh 


z 


Redd 


A r f 


-0.203 


0.057 


-0.201 


-0.033± 0.019 


0.394 ± 0.168 


radio- loud 


log A rf 


Redd 


2 


M BH 


-0.187 


0.218 


-0.224 


-0.225± 0.143 


0.950 ± 0.494 


radio-loud 


log z 


Redd 


A/bh 


A r f 


-0.153 


0.306 


-0.062 


-0.034± 0.116 


0.148 ± 0.028 


radio-loud 


log ibol 


z 


Redd 


A r £ 


-0.266 


0.297 


-0.269 


-0.058± 0.054 


2.794 ± 2.456 


radio-loud 


log L boi 


z 


M bh 


A r £ 


-0.384 


0.164 


-0.321 


-0.071± 0.055 


3.373 ± 2.536 


radio-loud 


log Redd 


z 


Z^bol 


A r £ 


-0.304 


0.300 


-0.274 


-0.045± 0.032 


0.079 ± 0.035 


radio-loud 


log Redd 


z 


M BH 


A r £ 


-0.369 


0.188 


-0.291 


-0.069± 0.047 


0.056 ± 0.055 


radio-loud 


log A/bh 


z 


Z^bol 


A r £ 


0.246 


0.386 


0.212 


0.040± 0.035 


-0.206 ± 0.309 


radio-loud 


log M B h 


z 


Redd 


A r £ 


-0.220 


0.432 


-0.180 


-0.052± 0.052 


0.600 ± 0.477 



Note. — Col. (1) The samples, namely the radio-quiet (radio-loud) subsamplc contains quasars with radio loudness smaller 
(larger) than 10. Col. (2) The parameter on which the dependence of variability is to be considered. Cols. (3-5) Parameters 
which are restricted in small ranges to exclude their influences on variability, denoted by xl, x2 and x3. Cols. (6-7) rl is the 
median value of Spearman rank correlation coefficients in all the qualified small data sets, which are obtained from binning 
in the restricted parameters, while p is the median value of the significance of its deviation from 0. Col. (8) rl is the median 
value of Pearson Product Correlation Coefficient, indicating the strength of linear correlations. Cols. (9-10) b is the median 
slope of the linear fitting for the dependence of the variability amplitude on the parameter denoted in Col. (2) and a is the 
median y axis intercept of the linear fit. 
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Fig. 1. — Normalized distribution of the variability amplitudes in the g', r' and i' band for 
the 416 radio- loud sources (red dashed lines) and 7242 radio-quiet sources (black solid lines). 
The median values are shown with the vertical lines. 
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Fig. 2. — Normalized distributions of the color differences between the brightest state and 
faintest state. There are 166 radio-loud (red dashed lines) and 2811 radio-quiet quasars 
(black solid lines) in the left panel, 158 radio-loud and 2553 radio-quiet quasars in the right 
panel, respectively. 
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Fig. 3. — log -Redd vs. z for the 7242 radio-quiet quasars. The bin size of -Redd and z 
are 0.5 dex and 0.3 respectively. Lines over-plotted divide quasars into 16x5 subsamples. 
The 49 grids filled with red points refer to the qualified subsamples and they are marked 
with algebraic numbers in the lower corner of each grid. They are further separated into 
sub-subsamples with log Mbh bin size of 0.5 dex. The algebraic numbers in the upper corner 
of the grids represent the indices of qualified sub-subsamples (with more than 10 quasars in 
each) for the study of the dependence of the variability amplitude on A r f. 
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Fig. 4. — Dependence of V on A r f in the g', r' and i' band. To demonstrate the relation, 
only that for the 25 largest qualified sub-subsamples obtained from Fig. [3] are shown. In 
each panel, the 10 red diamonds refer to the median values of V and A r f in the bins obtained 
from evenly dividing A r f of the sub-subsample. The legends in the upper left corner of each 
panel show the index of the sub-subsample obtained from Fig. [3] and the number of quasars 
it contains. rl_A r f, r2_A r f and 6_A r f are displayed in the upper right corner and ranges of 
restricted parameters are shown in the lower right corner. For the radio-loud quasars, refer 
to Fig. 21 in the Appendix. 
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Fig. 5. — Statistical results for the relationships between V and logA r f: p-\ r { vs rl_A r f(left 
panel), r2_X T{ (middle panel) and the linear fit results, |6_A r f |/|6 e rr-A r f | vs 6_A r f (right panel). 
Dots and solid lines refer to results of the radio-quiet subsamples, while asterisk signs and 
dashed lines refer to results of the radio-loud ones. 
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Fig. 6. — log .Redd vs. logMeH for the 7242 radio-quiet quasars. Their bin size are both 
0.5 dex. Among 5x6 subsamples divided by the overdrawn lines, 23 subsamples with more 
than 10 quasars in each (marked with algebraic numbers in the lower corner of each grid) are 
further divided into sub-subsamples with the A r f bin size of 400 A(see § El for more details). 
The algebraic numbers in the upper corner of the grids refer to the indices of the qualified 
sub-subsamples for the study of the dependence of V on 
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Fig. 7. — Dependence of V on z in the g', r' and z' band for the 25 largest qualified sub- 
subsamples extracted from Fig. |6j In each panel, the 10 diamonds refer to the median values 
of V and z in bins from evenly dividing z. The texts in each panel have the same meaning 
with Fig. HI For the radio-loud quasars, refer to Fig. 22 in the Appendix. 
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Fig. 8. — Figure similar to Fig. [51 but for the statistical results of the relationship between 
V and log z. 
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Fig. 9. — logMsH vs. z for the 7242 radio-quiet quasars. Their bin sizes are 0.5 dex and 
0.3 respectively. Among 6x16 subsamples divided by overdrawn lines, 52 subsamples with 
more than 10 quasars in each are qualified, represented by the grids filled with red points 
and marked with algebraic numbers. 
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Fig. 10. — logLboi vs. z for the 7242 radio-quiet quasars. Their bin sizes are 0.5 dex and 
0.3 respectively. Among 7x16 subsamples divided by overdrawn lines, 45 subsamples with 
more than 10 quasars in each are qualified, represented by the grids filled with red points 
and marked with algebraic numbers. 
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Fig. 11. — Dependence of V in the g', r' and z' filter band on Lboi for the 25 largest qualified 
subsamples obtained from the log -Redd - ^ space as shown in Fig. [3J In each panel, only V in 
the r' band for all the quasars are represented by grey points. The 10x3 colored diamonds 
(green, blue and red) show the median V values in the 3 bands (g', r' and i') in the bins 
obtained from evenly dividing log Lboi- The blue lines represent the linear fits for all the data 
in the r' band. The legends in the upper right corner of each panel list the rl_L b oi, r2_Lboi 
and fc_£boi for the study of the variability amplitudes in the r' band, while the texts in the 
lower right corner of each panel give the ranges of restricted parameters. For the radio-loud 
quasars, refer to Fig. 23 in the Appendix. 
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Fig. 12. — Figure similar to Fig. [TTJ but for the 25 largest qualified subsamples obtained 
from the log Mbh-z space as displayed in Fig. |9j Refer to Fig. 24 in the Appendix with 
respect to the radio-loud quasars. 
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Fig. 13. — Figure similar to Fig. [5l but for the statistical results of the relationship between 
V and logLboi- The magenta diamonds and triangles show respectively the results of the 
radio-quiet and radio-loud subsamples. These subsamples are only binned in redshift, thus 
without excluding the relations among the three intrinsic quasar parameters. The blue signs 
and lines represent results of the subsamples obtained from the logMBH-2 space and red 
ones refer to results of the subsamples binned in the log -Redd-- 2 space. 
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Fig. 14. — Dependence of V in the r 9 and i 9 filter band on -Redd for the 25 largest qualified 
subsamples obtained from the logLbor^ space as shown in Fig. [TD1 In each panel, only V in 
the r 9 band for all the quasars are represented by grey points. The 10x3 colored diamonds 
(green, blue and red) show the median V values in the 3 bands (</, r 9 and i 9 ) in the bins 
obtained from evenly dividing log -Redd- The blue lines represent the linear fits for all the 
data in the r 9 band. The legends in the upper right corner of each panel list the ?"1_Redd ; 
r2_i?EDD and 6_Redd for the study of the variability amplitudes in the r 9 band, while the 
texts in the lower right corner of each panel show the ranges of the restricted parameters. 
Refer to Fig. 25 in the Appendix for the radio-loud quasars. 
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Fig. 15. — Figure similar to Fig. dU but for the 25 largest qualified subsamples obtained 
from the logMBH-^ space as shown in Fig. [91 Refer to Fig. 26 in the Appendix for the 
radio-loud quasars. 
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Fig. 16. — Figure similar to Fig. [131 but for the statistical results of the relationships between 
V and log -Redd- The grey signs and lines represent results of the subsamples obtained from 
the log L ho \-z space and blue ones represent results of the subsamples obtained from the 
log Mbr-z space. 
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Fig. 17. — Dependence of V in the g', r' and filter band on Mbh for the 25 largest qualified 
subsamples obtained from the log -Redd - ^ space as shown in Fig. [3J In each panel, only V in 
the r' band for all the quasars are represented by grey points. The 10x3 colored diamonds 
(green, blue and red) show the median V values in the 3 bands (g', r' and i') in the bins 
obtained from evenly dividing log Mbh- The blue lines represent the linear fits for all the 
data in the r' band. The legends in the upper right corner of each panel list the rl_MeH, 
t2_Mbh and 6_Mbh for the study of the variability amplitudes in the r' band, while the texts 
in the lower right corner of each panel give the ranges of the restricted parameters. For the 
radio-loud quasars refer to Fig. 27 in the Appendix. 
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Fig. 18. — Figure similar to Fig. [T7J but for the qualified subsamples obtained from the 
log L^oi-z space as illustrated in Fig. [10J Refer to Fig. 28 in the Appendix for the radio-loud 
quasars. 
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Fig. 19. — Figure similar to Fig. [131 but for the statistical results of the relationship between 
V and logMeH- The red signs and lines represent the results from the subsamples obtained 
from the log -Redd - ^ space and grey ones represent the results from the subsamples obtained 
from the log L^ \-z space. 
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Fig. 20.— The upper left panel shows the model-predicted variability amplitude as a 
function of wavelength, here M B h = 10 8 ' 7 M and -Redd = 0.1259; the upper right panel 
shows the variability amplitude in r' band as a function of -Redd- The lower left panel shows 
the variability amplitude in the r' band as a function of Mbh; the lower right panel shows 
the variability amplitude in the r' band as a function of Lboi- The red, blue and black lines 
refer to the cases with different changes of accretion rate, namely 0.4m, 0.2m and 0.1m. The 
obtained slope b values for the simulated data during the case with the change of accretion 
rate at 0.2m are shown in the upper right corner of each panel. 
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A. Appendix 



For the radio-loud quasars, we use the same mechanism to investigate correlations be- 
tween V and quasar parameters. The detailed relations are shown here. 



-44- 



2_1 n_K„- -0.187 
n=15 <-2-X„= -0.153 
b_\„- -0.096 ± 0.094 

i I + + 

% ~tfc -|bo -eiSI.OO . 
+ -i.5 Tr™ a-i.o 

9.0 + M^Mj £ 9,5 


4_1 + 7fr_X„= -0.255 
n=14 J- r2-X„= -0.224 
Tjjl„- -0.225 ± 0.154 

■ 4 

T 1.00 <zS1.30 . 

-2.0 < R E „a-i.5 

9.0 < M„ (M^) £ 9,5 


5_1 rttA„- -0.086 
n=14 r2_X,i= -0.102 
-<rJi„- -0.093 ± 0.143 

++ + + + 

. 4 '- 00 <**'-30 . 
T -f.5 < r™ a-i.o 
-fl-.O < M„ (M^) S 9,5 


7_1 + rl-V= -0.272 
n=14 r2_X,i= -0.347 
b_X„- -0.316 ± 0.135 

v + 
A+ + 

- + f# 

+ -fl" 1.30 <zai.60 . 
-1.5< R™ a-i.o 
8.5 < M„ (Mo) i 9,0 


T 

7_2 _i_ -0.110 : 

n=14 r2_X„= -0.277 

buJk.-"PO»4 ± 0.161 

+ + + 

- H; : 

* I ++ " 

-f qt 1.30 <I*1.60 . 
-1.5< r eoo a-i.o 
9.0 < M„ (Mo) i 9.5 


a_i ri_v -0.122 

n=15 r2_A„ = -0.104 
b_A„- -0.184 ± 0.268 
+ 
+ + 

4r 4. 1 !-i.o < R EO n a-0.5 

8.5 < M m * 9.0 


10_1 rl_V -0.314 
n=12 r2_A„ = -0.320 
_|_ b_*,i- -0.236 ± 0.120 

+ 

: 

_L~^~i|- l.60<i*1.90 . 
M--q_TlL -i.o < r a-o.5 
~ 8.5 < M m (wQ) i 9,0 


10_2 r'Ai- -0.432 
n=14 r2_\,i= -0.409 
b_\,i- -0.256 ± 0.090 

+ 

T . _u_ 1.60 <zai.90 . 
T +"""-i.o < r™ a-o.5 
_|_ +.0 < M m (M^) i 9,5 


123456123456 



123456123456123456 

K, (10 3 A) 



Fig. 21. — Dependence of V on A r f in the g', r' and i' band for all the 8 radio- loud sub- 
subsamples within small ranges of -Redd? z and Mbh- In each panel, the red curve refers to 
the polynomial fitting with the function similar to Eq. (4), where X is A r f. rl_A r f, r2_A r f and 
6_A r f are shown in the upper right corner and ranges of the restricted parameters are also 
shown in the lower right corner. Refer to Fig. H] for the radio-quiet quasars. 
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Fig. 22. — Dependence of V on z in the g', r' and i' band and all the 18 radio-loud sub- 
subsamples within small ranges of -Redd, ^bh an d A r f. rl_z, r2_z and ranges of the restricted 
parameters are shown in the upper right corner of each panel. Refer to Fig. [7] for the radio- 
quiet quasars. 
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Fig. 23. — Dependence of V in the g', r' and filter band on Lboi in the qualified subsamples 
obtained from the logi?EDD--2 space, represented by green, blue and red plus signs respec- 
tively. The blue line refer to the linear fit for all the data in the r' band. The legends in 
the upper right corner of each panel refer to the rl_Lb i, r2_L^ i and 6_Lboi for the study of 
the variability amplitudes in the r' band and ranges of the restricted parameters. Refer to 
Fig. [TTJ for the radio-quiet quasars. 
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Fig. 24. — Figure similar to Fig. 23, but for the qualified subsamples obtained in 
logM B H-"Z space. Refer to Fig. [12] for the radio-quiet quasars. 
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Fig. 25. — Dependence of V in the g', r' and i' filter band on -Redd in the qualified subsamples 
obtained from the log Lt, Q \-z space, represented by green, blue and red plus signs respectively. 
The blue lines represent the linear fits for all the data in the r' band. The legends on the 
upper right of each panel list rl_R EDD , r2_R EDD and 6_Redd for the study of the variability 
amplitudes in the r' band and ranges of the restricted parameters. Refer to Fig. [TH for the 
radio-quiet quasars. 




Fig. 26. — Figure similar to Fig. 25, but for the qualified subsamples obtained in 
logM B H-"Z space. Refer to Fig. [15] for the radio-quiet quasars. 



- 50 - 



E 0.2 



rl_M w = -0.341 
4- r 2-M m = -0.228 
b_M~= -0.071 ± 0.092 
1 T 0.4Q < z S0.7O 
-1.5 < log R™*-I.O _ 

+ 


r1_M w = -0.222 
+ r2_M w = -0.247 
b_M w =j -0.067 ± 0.051 

~r o.70 < z a i.oo 

-1.5 < log f? t „*-1.0 _ 

+ + 


d-M»,= -0.251 
r2_M m = -0.014 
b_M w = -0.005 ± 0.092 
T 0.70 < z S 1 .00 
-1.0 < log R™*-0.5 _ 

4= + 


: fi_M w = -0.096 
4- r2-Mi,= -0.148 
b_M w = -0.057 t 0.098 

T i .oo 4-i a i .30 

-2.0 < log R t „5-1.5 _ 

"it* ' 


=1= r1_M*,= -0.128 
r2_M w = -0.215 
b_M w = -0.059 t 0.056 

4- i.oo < z a i.3o 

_j_ + - 1-jLiS log R e „*- 1 .0 . 


rl_M B1 = -0.188 
r2_M m = -0.358 
b_Mi w = -0.069 ± 0.042 
T 1.00 < z SI. 30 
+ 1-1.0 < log R £ „*-0.5 _ 
4+ 

£ 


f rl±M m = -0.135 
r2 r M B1 = -0.177 
b_M w = -03333 ± 0.033 
+ + T30 < z SI. EO 

_|_ ■qV*Hc log R tH1 *-l.O _ 


M_M„= -0.155 
r2_M B1 = -0.OB5 
b_M w = -0.034 ± 0.083 
+ 1.30 < z SI. 60 
_i_ +1-0 < log R £ „5-0.5 - 

^+++ - 




Vl_M w = -0.105 
r2_M w = -0.083 
b_M,„= -0.014 t 0.042 
+ 1.60 < z a 1.90 
-1.5 < log R £[o 5-1.0 _ 

++ + + 
+++ 

; 


r1_Ma,= -0.193 
r2_M w = -0.276 
b_M,„= -0.056 t 0.039 
1.60 < z SI. 90 


r1J« M » -0.709 
r2_M w = -0.659 
b_M w - -0.062 ± 0.021 
1.90 < z S2.20 
-0.5 < log Rojo* 0.0 _ 


8 9 10 8 9 10 8 9 10 8 9 10 



8 9 10 



log M BH (M ) 

Fig. 27. — Dependence of V in the g', r' and z' filter band on Mbh in the qualified subsamples 
obtained from the log R^dd-z space, represented by green, blue and red plus signs respec- 
tively. The blue lines show the linear fits for all the data in the r' band. The legends on the 
upper right of each panel list the rl_MBH> ?"2_Mbh and 6_Mbh for the study of variability 
amplitude in the r' band and ranges of the restricted parameters. Refer to Fig. [T7] for the 
radio-quiet quasars. 
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Fig. 28. — Figure similar to Fig. 27, but for the qualified subsamples obtained from the 
logL bol -2; space. Refer to Fig. [18] for the radio-quiet quasars. 



